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ABSTRACT 



We present the results of a wide-field survey using the 1.2-m Samuel Oschin Telescope 
at Palomar Observatory. This survey was designed to find the most distant members 
of the Kuiper belt and beyond. We searched ~12,000 deg 2 down to a mean limiting 
magnitude of 21.3 in R. A total number of 52 KBOs and Centaurs have been detected, 
25 of which were discovered in this survey. Except for the re-detection of Sedna, no 
additional Sedna-like bodies with perihelia greater than 45 AU were detected despite 
sensitivity out to distances of 1000 AU. We discuss the implications for a distant Sedna- 
like population beyond the Kuiper belt, focusing on the constraints we can place on the 
embedded stellar cluster environment the early Sun may be have been born in, where 
the location and distribution of Sedna-like orbits sculpted by multiple stellar encounters 
is indicative of the birth cluster size. We also report our observed latitude distribution 
and implications for the size of the plutino population. 

Subject headings: Kuiper belt: general- Oort Cloud 
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Introduction 



With the advent of wide- field CCD cameras in the past decad e, there has been an explosion in 



observational programs searching f or Kuiper belt objects fKBOs) (jJewitt fc Luu 



1996, 



1998; 



Elliot et al. 



Sheppard et al. 



200^ 



2000; 



Larsen et al. 



200 



Brown 



Larsen et al. 

r 



2001 



200 



Truiillo et a 



Kavelaars et al. 



2001 



1995; 



Jewitt et a 



Trujillo fe Brown 



2003 



2009). Now there are over 1000 



KBOs known, with about half having secure orbits. The majority of these surveys search for distant 
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solar system bodies using images taken on a single night over a span of a few hours, probing out to 
distances of ~100 AU. Most of these surveys have focused on observing within 10° of the ecliptic 
with the majority only imaging within just a few degrees. 



The discovery of Sedna (jBrown et al.l I2004J ) on a highly eccentric orbit far outside the Kuiper 
belt challenges our understanding of the solar system. With a perihelion of 76 AU, Sedna is well 
beyond the reach of the gas-gian ts and could not be scattered into its highly e ccentric orbit from 
interactions with Neptune alone ([Emel'vanenko et al.ll2003l ; lGomes et al.ll2005bl ). Sedna's aphelion 
at 1000 AU is too far from the edge of the solar system to feel the perturbing effects of passin g 
stars or galactic tides in the present-day solar neighborhood ([Duncan et al.lll987l ; lFernandeall997l ). 
Sedna is dynamically distinct from the rest of the Kuiper belt, and its unexpected discovery alludes 
to a population of icy bodies residing past the Kuiper belt with perihelia greater than 45 AU and 
semimajor axes greater than ~200 AU, beyond whic h Neptune is unable to raise the perihelia of 
scattered disk KBOs through resonant perturbations (JGomes et al.ll2005bl ). 



Sedna is the only body known to reside in this region. Sedna was found near perihelion at 
a distance of ~88 AU, at the motion limit and brightness limit of its discovery survey (Brown et 
al 2008). With one night imaging, pr evious KBOs surveys were likely insensitive to the objects 



in the Sedna region. To date, surveys (ILarsen et al. 



2007 



Brown 



2008 



Parker Sz Kavelaarsll2010l ) 



have been unsuccessful in finding additional Sedna-like bodies. In order to find the largest and 
brightest members of the Sedna population, we have been engaged in an observational campaign 
to survey the northern sky. We present the results of our search for distant solar system bodies 
covering ~12,000 deg 2 within 30° of the ecliptic. Rather than searching over a single night, we 
use a two-night baseline to distinguish the extremely slow motions of these distant bodies from 
background stars. We are sensitive to motions out to a distance of ~1000 AU (~0.2 "hr _1 ). 

In this paper, we discuss the implications for a distant Sedna-like populatio n beyond the 



Kuip er belt and provide constraints on the cluster birth Sedna formation scenario (JBrasser et al. 



2006). The survey was specifically designed to find the select brightest members of a distant Sedna 
population but was also sensitive to the dynamically excited off ecliptic populations of the Kuiper 
belt including the hot classicals, resonant, scattered disk, and detached Kuiper belt populations. 
We present our observed latitude distribution and implications for the plutino population. 



Observations 



Observations were taken nightly using the robotic 1.2 m Samuel Oschin Telescope located 
at Palomar Observatory and the QUEST large-area CC D camera. The QU EST camera has an 
effective field of view of 8.3 deg 2 with a pixel scale of 0.87" ( Baltay et al.ll2007l ). The 161-megapixel 
camera is arranged in four columns or "fingers" along the east-west direction each equipped with 28 
2400x600 CCDs in the north-south direction (see Figure H]). The gap between CCDs in the north- 
south direction is ~1.2' and the spacing between adjacent fingers along the east-west direction is 



~25'. The four fingers are labeled (A-D) and the CCDs are numbered sequentially (1-28) from 
North to South. We will refer to the CCDs by finger and position along the finger (i.e.,C14, D28). 

Observations were taken from 2007 May 8 - 2008 September 27. We have surveyed in total 
11,786 deg 2 within ±30° of the ecliptic to a mean depth of R magnitude 21.3. Our sky coverage 
is shown in Figure EJ Field centers are compiled in Table [Do A forest fire on Palomar Mountain 
prevented observations in 2007 September and camera malfunctions ceased operations from 2008 
February -2008 May leading to gaps in longitudinal coverage. After 2008 May normal observations 
resumed until the QUEST camera ceased operations on the Oschin telescope at the end of 2008 
September. 

Target fields were observed over a two-night baseline in order to search for solar system objects 
out to distances of ~1000 AU (moving at speeds as low as 0.2 "hr _1 ). All exposures were taken 
through the broadband re d RG610 fi l ter ( IIIaF filter from the POSS-II survey) with a wavelength 



range of A=610— 690nm (jReid et al.l 119911 ). For each field, a pair of 240s exposures was taken 
separated by ~1 hour on each of the two nights. The second night of observations was typically 
the next day or at most four nights later. Observations were in varying photometric conditions 
and lunations. To check the photometric quality of each nightly pair of observations, magnitudes 
of the detected sources from both images were histogram binned with a bin size of 0.2 mag, and 
the peak value of the histogram was selected as an indicator of image depth. If the median value of 
the five CCDs best CCDs (B11,C19,D09,D12,D13) was less than 20.4 mag (19.0 mag for crowded 
fields with greater than 4000 detected sources) than the observation was rejected as poor quality, 
and the target field was rescheduled for new observations the next night. If a target field cannot 
be successfully imaged within four nights of the first pair of observations, the field was reset and 
scheduled for another two nights of observations. 

All target fields were observed within 42° of opposition, and to avoid high star densities, fields 
less than 15° from the galactic plane were avoided. The camera RA CCD gap was covered by adja- 
cent pointings, but the ~1.2' declination gap remains mostly uncovered in our survey observations. 
When all opposition fields within ±30° of the ecliptic for a month's lunation were completed, over- 
lap pointings were then targeted to reduce holes in our sky coverage due to the camera's declination 
gap and defective CCDs. From the beginning of the survey to 2007 November 12, instead of per- 
forming overlapping coverage, fields with ecliptic latitudes greater than 30° were instead targeted 
once all available opposition fields within 30° of the ecliptic were completed. 



the full version of Table 1 will be available in the online journal version 



3. Data Analysis and Object Detection 

3.1. Moving Object Detection 

Observations were processed nightly though an automated reduction pipeline using the In- 
teractive Data Language (IDL) Software package. Each CCD on the detector was reduced and 
searched for moving objects independently from the other CCDs on the mosaic. All images were 
bias-subtracted and flat-field corrected. A row- by-row median of the overscan region was used for 
the bias subtraction. A master flat-field image for each CCD was constructed from a 3-c clipped 
median of the night's science images. Some of the camera CCDs had a significant fraction of hot 
or defective pixels. These pixels were identified as those where the flat field image value deviated 
by more than 0.7% from the value of the 3x3 median boxcar filtered flat field image. To mask the 
effects of these hot pixels, those regions of the science image were replaced by the median value 
of a 3x3 pixel box centered on the bad pixel. If hot/bad pixels constituted more than 20% of the 
image than a 3x3 median boxcar smoothed image was substituted for analysis. 



SExtractor (jBertin fc Arnoutslll996l ) was run on each image to compile a list of sources. SEx- 
tractor was tuned such that source detection constituted 4 or more contiguous pixels (DETECT_MINAREA 
parameter) above the detection threshold (DETECT_THRESH parameter) of 1.2 a above the sky 
background. Chips C14, C04, and D26 had significant image defects and higher SExtractor detec- 
tion thresholds were used for these three CCDs (DETECT_MINAREA =5 and DETECT.THRESH 
=2.3) SExtractor performed circular-aperture photometry using a 5-pixel radius aperture. Each 
source was characterized by its position, flux, and shape. The best of the four images was selected 
as the master template wh ose astrometric solution was found by matching image stars to the USNO 



A2.0 catalog (|Monetlll998l ). The other three images were then aligned relative to the stars in the 
template image. Even if the absolute astrometry failed, the relative astronomy between the images 
was still sufficient to search for distant solar system bodies. The median absolute astrometric error 
for the entire survey was 0.4". The median relative astrometric error between survey images was 
0.076". 

Once astrometric solutions had been found, the images were searched for moving objects. 
Because our observations were taken at or near opposition, slow-moving solar system objects were 
identified by their retrograde motion due to the parallax caused by the Earth's orbital motion. 
Distant planetesimals may move too slowly to show apparent motion over the nightly one-hour 
baselines and appear stationary on individual nights. To ensure the detection of objects out to 
distances of ~1000 AU, we only required motion to be identified over the two-night baseline. The 
detection catalogs from all four images were compared to identify and eliminate the stationary 
sources in each image. Sources on one image that had a counterpart within a 4" radius on either 
of the second night's observations were removed as background stars. To further cull the object 
lists of stars that were above the SExtractor thresholds on one night but below the detection limit 
on the other, we generated SExtractor source catalogs with more sensitive detection parameters 
(DETECT_MINAREA = 3 and DETECT.THRESH = 1.1), and compared these deep catalogs to 



our detection lists. Image sources from one night that appeared on the other night's deep detection 
catalogs were deemed stationary and rejected as well. Saturated stars and extended sources whose 
peak flux was more than 3 pixels from the source center measured by SExtractor were also removed 
from the object catalogs. 

Potential moving candidates were identified from the remaining unmatched sources. The 
nightly images were searched for moving object pairs with motions less than 14.4 "hr _1 , the velocity 
of bodies at distances of 10 AU. Moving object pairs from the first night and pairs from the second 
night separated by more than 4.38" with retrograde motion consistent with opposition were linked. 
To eliminate stationary image sources that had been linked between the two nights, candidates 
with average nightly magnitudes differing by more than one magnitude were eliminated. Remain- 
ing candidates whose nightly motions differ by less than twice the first nights on sky velocity were 
kept to create the lis t of moving object candidates. Candidates were filtered via the orbit-fitting 



package described in iBernstein fc Khushalanil (|2000l ) . Those candidates with successful orbit fits 
which produced a x 2 l ess than 25 and barycentric distance between 15 and 1000 AU were identified 
as moving objects and added to the final list of candidates to be screened by eye. 100x100 pixel 
subimages for each of the final moving object candidates were created from the discovery images. 
These snapshots were aligned and blinked by eye. A total of 39,110 candidates (~200 a night) 
were visually inspected. Typical false positives included diffraction spikes, faint background stars, 
blended sources, and CCD imperfections. 



3.2. Recovery Observations 

At discovery, heliocentric distance and inclination can be identified from the parallax effect 
due to the Earth's motion, but other orbital parameters remain unconstrained. With only a two- 
night discovery arc, a distant Sedna-like body cannot be distinguished from a typical scattered disk 
Kuiper belt object near aphelion. Even with follow-up observations a month after discovery, both 
families of orbits provide reasonable astrometric fits to the observations. The two orbital solutions 
diverge sufficiently a year after discovery, and a secure dynamical identification can only be made 
after these additional observations. 

Recovery observations of new discoveries were taken at the Palomar 60-inch telescope, the 
Palomar 200-inch telescope, the 0.9-m telescope operated by the SMARTS consortium at Cerro 
Tololo Inter- American Observatory , the 42-inch John S. Hall Telescope located at Lowell Observa- 
tory, the 2.66-m Nordic Optical Telescope located at el Roque de los Muchachos Observatory, and 
then 8.2-m Subaru Telescope on Mauna Kea. Of our detected KBOs, 96% have multi-opposition 
observations. All but two discoveries classified as KBOs by the Minor Planet Center (2007 JF45 
and 2007 PS45) were recovered during the survey. The two unrecovered objects were discovered 
during reprocessing of the data with more sensitive SExtractor source detection parameters and 
were discovered after they were no longer observable. Observations taken near 40° from opposition, 
contained contamination from asteroids near their stationary points that appeared to be moving at 



rates similar to distant KBOs. Some were identified with subsequent observations that confirmed 
these objects were on orbits with semimajor axes less than 5 AU. All other objects not successfully 
recovered have either been linked with other asteroid observations or have been classified on orbits 
well short of the Kuiper belt by the Minor Planet Center (MPC) database^. 



3.3. Calibration and Efficiency 

3.3.1. Limiting Magnitude 

The survey observations were taken during a wide variety of photometric, seeing, and weather 
conditions. Each CCD frame was independently photometrically calibrated. A photomet ric zero 



point offset to our instrumental magnitudes was derived relative to the USNO A2.0 catalog ((Monet 



19981 ) red magnitude. The photometric uncertainty of th e USNO cata log is non-negligible. For 



magnitudes greater than 17, the uncertainty is 0.3 mag ([Monet] 119981 ). We likely have several 
tenths of magnitude uncertainty in our discovery magnitudes. We have not precisely calibrated 
the survey depth with calibration observations. Limiting magnitudes were computed based on the 
USNO catalog. We found that the faintest magnitude with a 5a (10a for C2 A19, C14, C04, D26; 
CCDs with larger numbers of hot pixels), uncertainty as reported by SExtractor represented an 
accurate measure of the source detection limit of our images, and we used these values in the work 
presented in this paper. The limiting magnitude for each nightly pair of field observations was 
taken as the depth of the shallower of the two images. The mean limiting magnitude of the survey 
based on the USNO catalog is 21.3 in R. 



3.3.2. Survey Efficiency 

Because our survey has covered a wide swath of sky detecting multiple previously known 
KBOs, we have an alternative method of determining the limiting magnitude of our survey. Of 
our detections, 27 are previously discovered KBOs and Centaurs in the MPC database. The 
absolute magnitudes recorded in the MPC are based upon the apparent magnitudes measured from 
the discovery or follow-up observations, like our survey, which are often taken in n o n-stan dard 



filters and observed without precise photometric calibrations. iRomanishin &: Teglerl (|2005l ) find 
the absolute magnitudes recorded in the MPC are systematically 0.3 mag brighter than those 
magnitudes accurately measured for their sample of 90 KBOs and Centaurs. We can still use the 
known population of bright KBOs to estimate a crude efficiency for the survey. We obtained the 
positions and visual apparent magnitudes computed by JPL Horizonajj for known KBOs. As of 2010 



2 http : //www . cf a . harvard . edu/iau/Ephemer ides/Distant /index . html 
3 http : //ssd . jpl . nasa . gov/horizons . cgi 



January 20, there were 64 previously known multiopposition KBOs with visual magnitudes brighter 
than 22nd magnitude (excluding discoveries found in this survey and objects with a < 30) with 
predicted positions located on our survey images that could have been detected by our detection 
pipeline. We only considered KBOs positioned on the same CCD for all 4 field observations, not 
accounting for masked regions of the CCDs. Masked bad pixel regions account for ~8% of the 
QUEST camera's observable area, but likely the loss due to bad pixels is much smaller than this 
value. A KBO positioned on a bad pixel may not necessarily be lost, Sextractor interpolates values 
for masked pixels from neighboring good pixels before source detection. For every object not 
detected in the survey, we examined the images to determine if a moving source was visible. No 
known KBO was missed during the visual inspection of moving object candidates. The majority of 
the missed KBOs were not found because the KBO's psf overlapped with a neighboring star and 
was missed by SExtractor, the KBO was on a bad or masked off region of the CCD, image quality 
was bad due to poor telescope tracking, or the KBO was too faint to be detected and no visible 
moving source was identifiable. 

We define the survey efficiency function as: 

1 - tanh J™' 171 * ) ) (1) 



"2 V V 9 

where e is the efficiency with which KBOs of magnitude m are detected in our survey, e max is 
the maximum efficiency, m* is the magnitude at which £m ^ x , and g is the half width. We fit 
for the efficiency by computing the cumulative distribution for all known KBOs scaling for the 
probability of detection and compare to the observed cumulative distribution. To find the optimal 
parameters, we minimize the x 2 between the observed and calculated cumulative distributions. We 
find e max =0.6Q, m*=21.5, g=0.05. The efficiency drops by 50% at 21.5 V mag, consistent with 
our median image limiting magnitude. Figure [3] plots the best-fit efficiency function and plots the 
binned detection efficiency for the known sources located on all 4 field images in 0.5 mag bins. 
We estimate the uncertainty in our survey efficiency using the number of known KBOs found with 
magnitudes less than or equal to 21st magnitude, well before the drop off in the best-fit efficiency 
function. We found 13 of 19 known KBOs brighter than or equal to 21st magnitude, giving an 
efficiency of 68%, consistent with our best-fit efficiency function, and assuming Poisson counting 
statistics, the 1-cr confidence level ranges from 51 



3.3.3. Geometric Losses 

The gap between the QUEST camera's CCDs in the north-south direction is ~1.2'. Along the 
east- west direction the separation between CCDs is ~25'. At the distances our survey is sensitive 
to, a KBO located in the declination gap would remain in the gap between CCDs over the two-night 
baseline. This was the case for Eris, and Eris was not detected in our survey. For some areas of the 
sky we do have overlap pointings to try and cover the declination gap but only after all opposition 



target fields were observed. The losses due to the CCD gaps is accounted for in our sky coverage 
estimates, but we do not include the loss due to masked regions. 

KBOs that moved off the edge of the CCD into the CCD gaps were missed by our automated 
detection pipeline. Non-functioning CCDs and longitudinal losses are accounted for in our lati- 
tudinal sky coverage estimates, but to measure our geometric losses from those KBOs moving off 
the CCDs or lost in the CCD gaps, we generated ~10 6 random circular orbits assuming a uniform 
inclination distribution (0-180°) for a range of semimajor axes. Neglecting the effect of masked 
CCD regions, we calculated the fraction of simulated KBOs positioned on all four survey images as 
a function of ecliptic latitude. Figure H] compares our survey sky coverage to the fractional coverage 
of the simulated circular orbits at 30, 50, and 100 AU. The greatest losses occur at the ecliptic, 
and we find this effect is at most ~10%. Closer orbits are moving at faster on-sky velocities and 
are more likely to move off the CCD over the two night-base line than objects at further distances, 
but we find the difference in losses by objects at 30 and at 100 AU is small, and that all objects in 
the Kuiper belt have similar geometric losses in our survey. 



3.3.4- Pipeline Detection Efficiency of Sedna-like Bodies 

Any comparison of the Sedna population requires that we also understand whether these bod 
i es would be detected in our survey. Many of the mechanisms proposed for the formation of Sedna 



(Ke nvon &: Bromley 



2004 



Morbidelli fc Levison 



2004 



Brasser et al. 



2006 



20071 : iKaib fc Quinn 



2008) produce many highly eccentric and even retrograde orbits. To test whether Sedna-like orbits 
would pass through our orbit-fitting filter, we created artificial orbits with a uniform semimajor 
axis ranging from 100-1100 AU and uniform eccentricity and inclination distribution including 
retrograde orbits. For those 781,763 artificial orbits whose positions land on our images, have 
barycentric distances less than 1000 AU and have perihelia greater than 50 AU, we add absolute 
and relative positional offsets characteristic of the survey's astrometric errors. All four images of 
a field observation have the same absolute astrometric error but random relative positional errors. 
We add normally distributed random absolute and relative astrometric errors using the three-sigma 
clipped median and standard deviation of the survey astr ometric uncertainties. As shown in Fig [5j 
the efficiency is the fraction of synthetic orbits fit with the iBernstein fc Khushalanil ()2000j) software 
that pass our selection criteria in each semimajor axis bin compared to the number of objects in 
the 100 AU bin. 5% of the simulated orbits would not have made it through to visual inspection 
with the majority of failures due to the best-fit orbit placing the object on an asteroid-like orbit. 
We are confident that Sedna-like bodies present in our images detected by SExtractor would be 
identified by our automated detection scheme. 



Detections 



A total of 52 KBOs and Centaurs have been detected of which 25 are new discoveries from 
this survey. 50 of our discovered objects have multiopposition orbits. Table [2] lists the orbital 
information for objects detected in the survey. The orbital and radial distribution is plotted in 
Figure [6] and Figure [7] respectively. With the detection of no cold classical belt objects; our overall 
survey probes the orbital properties of the hot classical, scattered disk, detached, and resonant 
populations. The survey was specifically designed to probe the Sedna region, but except for Sedna, 
no additional objects with perihelion greater than 45 AU were detected despite sensitivity out to 
distances of 1000 AU. 

We do detect several Centaurs with semimajor axes less than 30 AU in our survey, but to 
constrain the number of false detections by inner solar system objects we placed a minimum distance 
threshold in our moving object detection scheme. C andidates with barycentr i c dist ances less than 
15 AU as calculated from initial orbits fit by the iBernstein fc Khushalanil (|2000l ) method were 
ignored. Our survey is limited to detecting only the most distant of the Centaurs, and we therefore 
will not address the Centaur population in this paper. 



5. Sedna Population 

With a perihelion of 76 AU and an aphelion of ~1000 AU Sedna is dynamically distinct 
from the rest of the Kuiper Belt. Its extreme orbit suggests the presence of a population of icy 
bodies residing past the Kuiper belt. The study of this Sedna population provides a unique new 
window into the history of the early solar system. Some other mechanism no longer active in 
the solar system today is required to emplace Sedna on its highly eccentric orbit. Several possible 
scenari os have been offered to explain Sedna's extreme orbit, including interactions with planet-size d 
bodies dciadman fc ChanlhoOfil: lOomes et all20od : lLvkawka fc Mukailboosl : IComes fc Soareslhoioh . 



stellar encounters (IMorbidel 



Morbidelli fc Levisonl booi 



Kenvon &: Bromlevl 



2004 



companion (JMatese et al. 



i & Levisor 



nl l2004h. 



Brasser et al] 



2006, 



multiple stellar fly-bys in a stellar birth cluster 



2007 



Kaib h Quinnl |2008| ) , interstellar capture 



Morb idelli fe Levisonll2004l ). and perturbations from a wide-binary solar 



20051 ). Each of the various Sedna formation models leave a distinctive 



imprint on the members of this class of distant objects and has profound consequences for our 
understanding of the solar system. These planetesimals in the Sedna region are dynamically frozen 
and the relics of their formation process. The orbital distribution and number density of Sedna-like 
bodies will distinguish between the formation scenarios. 



In ISchwamb et al.l ([2009), before recovery observations were complete, we compared the ex- 
pected number of detections from a theoretical population on orbits with the same semimajor axis 
and eccentricity as Sedna to our survey results, the redetection of Sedna. Our best-fit value gives 40 
bodies residing on Sedna's orbit that are brighter than or equal to Sedna. At the one-u confidence 
level we ruled out a population larger than 92 and smaller than 15 Sedna-sized or bigger objects on 



10 



orbits similar to Sedna's. Our previous work had been limited to examining a model population of 
bodies residing specifically on Sedna's orbit. Any realistic Sedna population likely occupies a much 
larger region of orbital space, possibly including objects with sufficiently high perihelia that they 
would never or rarely become bright enough to see. With secure orbital classifications for survey 
objects, we can now test more sophisticated orbital distributions. 



5.1. Constraints on a Cluster Birth 



No new Sedna-like bodies with perihelia beyond 45 AU were found in the survey despite a 
sensitivity out to distances of ~1000 AU. Although, we cannot differentiate between the Sedna 
origin scenarios with a single detection, we can place constraints on the cluster birth model where 
the location and distribution of Sedna-like orbits is indicative of the Suns birth cluster size. Most 



stars are born in dense gas-rich embedded clusters (JLada et al. 



2003 



Lada fe Lada 



2003 



1991 



Carpenterll2000l ; iPorras et al. 



Allen et al.ll2007l ). and it is likely that the Sun spent several million years 



in such an environment. The presence of short-lived radioactive nuclides in primitive meteorites, 
may provide circumstantial evidence that the Sun was in relatively close proximity to a supern ovae 
early on in the solar system's formation, (jChaussidon &; Gounell e 2007 B rennecka et al.ll2009i and 
references therein) and therefore in a much denser environment than the present-day solar neigh- 
borhood. In the dense stellar nursery, encounters between nearby solar nei ghbors and the Sun 
woul d occur at a much higher f r equency than in the present solar environ ment ([Adams &: Laughlin 
200 ll : lLaughlin fc Adams! 1 19981 ; IProszkow fc Adamsl 120091 ; Udamslboiol ). Close fiy-bys of passing 
stars would perturb objects in the Sun's planetesima l disk onto highly ecc entric Sedna-like orbits 
(|Morbidelli & LevisonlEool : brasser et aObood . I2OO7I : kaib & Quirmlhoosl ). 



Brasser et al.1 (120061 ) successfully produce objects on orbits similar to Sedna's in simulations 



of embedded cluster environments. The gravitational effects of both stars and gas in the cluster 
are included in their integrations. If the mean density of the material the Sun encounters while 
residing in the embedded cluster was ~10 3 Mq/pc 3 (central cluster densities of 10 4 Mq/pc 3 ) or 
denser, Sedna's orbit is rec reated and a d i stribu tion of Sedna-like bodies with semimajor axes less 
than 10,000 AU is formed. iBrasser et al.l (|2006l ) find that the central density of the stellar cluster 
(directly correlated to the amount of material the Sun encounters in the cluster) determines the 
orbital distribution of Sedna-like bodies generated. The denser the cluster environment, the smaller 
semimajor ax i s at w hich the Sedna population begins. For this paper, we focus specifically on the 
Brasser et al.1 (120061 ) results for the 10 4 ,10 5 , and 10 6 Mq/pc 3 embedded cluster integrations (10 3 
M /pc 3 did not produce Sedna). We refer the read er to their pa p er for details of the orbital 



integrations and the review of embedded clusters by lLada fc Ladal (120031 ) . Figure [8] shows the 
orbital distributions from the embedded cluster numerical simulations used in this work. 

Our survey observations probe the Sedna po p ulatio n today after 4.5 Gyrs of evolution. The 
distribution of orbits presented by IBrasser et al.l (J2006J) is what remains after 3 Myr when the 
integrations end and the Sun is expected to have left the birth cluster. Once the Sun exits the 
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cluster, the Galaxy becomes the dominant gravitational potential. The gravitational perturbations 
from galactic tides over the age of the solar system have not been accounted for in the lBrasser et al 



(J2006I ) integrations. Sedna's orbit is protected from the effects of passing stars and galactic tides 
in the current solar environment, but objects with h i gher semimajor a xe s than Sedna may be 
perturbed onto comet-like orbits (JDuncan et al.1 119871 ; iFernandea 119971 ). iKaib fc Quinnl (|2009l ) 
examined the production of long period comets in the Sedna region and find that the production 
efficiency drops significantly for bodies with a < 3000 AU compared to those with larger semimajor 
axes. Therefore we expect that objects emplaced onto Sedna-like orbits with semimajor axes less 
than 3000 AU should remain to the present day, and we do not include any orbits from the cluster 
simulations with a > 3000 AU in comparisons to our observations. 



The Nice model (jTsiganis et al. 



2005 



Gomes et al.ll2005al ) predicts that the giant planets were 



in a more compact configuration than in the present-day solar system. The orbits of the giant 
planets went unstable approximately 1 Gyr after the formation of the solar system causing the 
migration of the giant planets and scattering of planetesimal disk. Jupiter migrates inward, and 
the remaining giant planets move ou tward with Neptune migrating outward t o 30 AU. The oldest 
embedded clusters are ~5 Myrs old (JLeisawitz et al.l 1 19891 ; lLada k, Ladal 120031 ) , Neptune migrates 
well a fter the Sun has left the birth cluster and the emplacement of the Sedna population. iBrasser 
(|2008l ) confirms this scenario can create a Sedna population and generate an Oort cloud population 
within the current estimates of the mass of the Oort cloud. Neptune's orbit became eccentric during 
migration and was later circularized via scattering of planetoids in the Kuiper belt region. Current 
estimates have Nept une's eccentricity as high as ~0.3 corresponding to an aphelion of ~39 AU at 
the end of migration (jLevison et al.ll2008l ). The sculpt ing of the Sedna pop ulation due to Neptune's 
migration outward has not been accounted for in the lBrasser et al.l (120061 ) simulation results. The 
cluster models do create orbits with perihelia in the range of 30-50 AU, which may not exist in the 
current solar system due to Neptune ejecting these Sednas or scattering them onto KBO-like orbits 
during its eccentric phase. We chose a conservative minimum perihelia threshold of 50 AU (which 
would require Neptune to have an eccentricity of ~0.7 to reach 50 AU at aphelion) to compare the 
cluster distributions to our survey results. 



5.2. Survey Simulator 

We developed a survey simulator to compare the expected number of detections from the the- 
oretical cluster Sedna populations to our survey results. The simulator draws synthetic objects 
from a model orbital and absolute magnitude distribution and for every image computes the po- 
sitions and brightnesses of these objects on the sky. For all three cluster environments, we model 
a population of 3,000,000 bodies on cluster-created orbits randomly drawing the semimajor axis, 
eccentricit y, and inclina t ion fo r each particle from those produced in the cluster numerical inte- 
grations. IBrasser et al.l (|2006l ) obtain a value of ~2 Gyr for the precession frequency of Sedna 
and other Sedna-like objects, therefore we assume that the orbits have been randomized due to 
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planetary effects, and randomize over all other orbital angles. The positions of the artificial ob- 
jects are computed for the survey period; those synthetic cluster objects that land on our images 
are identified neglecting the effects of masked regions of the CCDs. For each of the three clus- 
ter environments, approximately a third of the synthetic cluster-created orbits are located on our 
images. 

A brightness distribution is then applied to the synthetic population. Due to the large uncer- 
tainties in the albedo distribution of such a distant population, we assign absolute magnitudes to 
our synthetic bodies instead of diameters. We assume a single power-law brightness distribution 
similar to the Kuiper belt where the number of objects brighter than a given absolute magnitude, 
Umax, is described by: 

N(tf < H max ) = iV H <i.6lO a(H — " L6) (2) 

The brightness distribution is scaled to Nh<\.q, the number of bodies with an absolute magnitude 
brig hter than or equal to Sed na (H=1.6). A typical value of a measured for the asteroid belt is 



0.3 (jJedicke &: Metcalfe! 1 19981 ). The best-fit single power law for the hot (inclinations > 5°) and 
cold (inclinations < 5°) Kuiper belt populations are a=0.35 and a=0.82 respectively, measured 
by Fraser et al (2010), but it is unclear if the Sedna population should have an a value similar to 
the Ku iper belt. The Sedna population may have very different surface characteristics than typical 



KBOs. iBarucci et all (120051 ') find methane and a tentative detection of nitrogen on Sedna's surface. 
Schaller fc Brownl (|2007l ) model of volatile loss on KBO surfaces predicts that moderate-sized Sedna- 
like bodies on high perihelia orbits should retain methane and nitrogen ices on their surfaces. Most 
KBOs on the other hand, are either too small or too hot to hold on to their primordial abundance 
of volatiles. Distant Sednas never sublimate a significant amount of ices to renew their surfaces 
in a frost /thaw cycle. Instead the surfaces of the Sedna population would be subject to constant 
photoprocessing of methane by solar irradiation steadily darkening their surfaces. Sedna is one of 
the reddest KBOs with a V-R=0.78 w i th thermal measureme nts constraining Sedna's V albedo to 



be between 0.16 and 0.30 (JBrownll2008l ; IStansberrv et al.ll2008l ). We choose to explore the extremes 



of the brightness distribution and model the likely range of power-law distributions for a ranging 
from 0.2-0.8 including the best-fit value for the hot and cold KBOs measured by Fraser et al (2010). 

For a given value of a and Nh<i.g absolute magnitudes are randomly assigned to our simulated 
Sednas. A single instance of the brightness distribution can be thought of as a separate survey. 
Those synthetic objects that lie within our sky coverage with an apparent magnitude above both 
nights' limiting magnitudes (as determined in Section 13.3. ip are deemed valid survey detections. 
We assume a 100% efficiency out to the limiting magnitude where then the efficiency immediately 
drops to zero. We require that the object must be located on all 4 field images to be considered 
"discovered" in the simulated survey, and we do not require the object have Sedna's perihelia of 
76 AU. Bodies with H < 4.3 residing at 50 AU would be visible within our survey, and an object 
of Sedna's size and albedo would have been detected up to a distance of ~93 AU. Objects have 
multiple detection opportunities due to repeat sky coverage over subsequent years and overlapping 
fields. We do not count duplicate detections in our tallies. 
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5.3. Could Sedna have been formed in a cluster environment? 

We did not find any distant objects with perihelia greater than 45 AU with the exception of 
Sedna. To determine whether the orbital distributions produced in the various cluster environments 
are consistent with our redetection of Sedna, we must compare the orbital distributions of single 
detections produced by the survey simulator to Sedna. We employ our 3,000,000 synthetic Sedna 
population for each cluster environment to generate single detections. For each given value of a, 
absolute magnitudes are randomly assigned to our simulated Sednas for the range of possible values 
of Nn<i.e to create 10,000 single detections. 

Each simulated detection is characterized by a semimajor axis, eccentricity, and inclination ( 
a,e,i). We test a,e,i because these three parameters are directly effected by the impulses from 
the stellar encounters and gravitational effects from the embedded gas and stars, and these are 
the most independent set of orbital parameters. We choose to exclude the H distribution in our 
analysis because of the uncertainty of our limiting magnitudes. To determine whether Sedna 
and the cluster produced single detections could be drawn from the same parent population for 
varying slopes and scaling of the brightness distrib ution, we emplo y a variant of a 3 -dime nsional 



( 3-D)Kolmogorov-Smirnov (KS) test adapted from iPeacockl (|1983l ) and lPress et al.l (|1992l ) which 
simultaneously compares the a,e,i orbital distributions to Sedna. The fraction of data points in 
each of the 8 quadrants in a, e, i space, where the origin is defined by Sedna's orbital parameters 
(a=519 AU , e=0.853, i=11.9 deg), is computed. In order to determine if Sedna's orbit is extreme 
compared to the cluster produced detections, the D statistic in this case is defined as the difference 
of the maximum and minimum fraction calculated. The significance of the computed D statistic is 
found by performing our 3-D KS test again, selecting each of the 10,000 simulated single detections 
as the new origin, counting the fraction where the computed D statistic was higher than the D 
statistic for Sedna's orbit. We reject the cluster-produced population if the 3-D KS test does reject 
at a 95% or greater significance the null hypothesis, that the simulated survey single detections 
and our sole detection of Sedna are drawn from the same distribution. 

We performed the 3-D KS test for all ranges of Njjkis that produced single detections and 
possible values for a (0.2-0.82) for all three cluster environments. The orbital distribution of single 
detections produced at smaller Nh<i&, is different from those at large Nh<i.6i an d the entire 
range of possible values Nh<is must be tested. At small values of Nh<i.g there are fewer bright 
H objects available to fill detectable orbits, biasing the single detections to slightly lower perihelia 
orbits than for larger values of Nh<i.q where there is an ample supply of bright bodies to fill 
detectable orbits. We find that the 3-D KS test confidence levels calculated for the 10 6 and 10 5 
M Q /pc 3 cluster distribution for varying values of a are independent of Nh<is- For the 10 4 cluster 
and any value of a, Nh<is = ^- nas the highest probability of rejection and then decreases to a flat 
value as Nh<i.g increases. For the 10 4 cluster, Nh<i.6=1 represents an upper limit on the rejection 
confidence level of the orbital distribution. Therefore we report the confidence level calculated for 
each cluster distribution and brightness distribution for values of Nh<is = ^ in Table O For the 
two densest cluster environments 10 6 and 10 5 M /pc 3 producing Sedna as the sole detection is 
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an extremely low probability event. The bulk of the 10 6 and 10 5 M Q /pc 3 cluster-created single 
detections had orbits with semimajor axes less than Sedna's. The simulations produce many more 
objects with lower perihelia than Sedna that should have been found but were not detected in our 
survey. We can rule out 10 6 and the 10 5 M /pc 3 cluster population at confidence levels greater 
than 95% for all ranges of a and possible values of Njj<i.6- Therefore we reject the 10 6 and 10 5 
M /pc 3 clusters as the source of the Sedna population. We cannot reject the 10 M /pc 3 cluster 
environments to a confidence level greater than 60% for all combinations of a and Nhkis tested; 
we are unable to rule this population out with statistical significance. The 10 M Q /pc 3 orbital 
distribution is consistent with our redetection of Sedna. These results assumed that every object 
that lands on a CCD brighter than the image limiting magnitude would be detected. Our detection 
efficiency is not 100%, but including a flat detection efficiency curve that drops to zero at the 
image limiting magnitude does not change the results presented. Including a detection efficiency 
produces the same types of orbits for single detections, just the absolute number of single detections 
decreases. Since we are only looking at single detections, the 3D KS test results are the same for 
any efficiency value. 

If Sedna's orbit is the result of multiple stellar encounters when the nascent Sun resided in an 
embedded cluster, our work rules out central densities for th e cluster greater tha n or equal to 10 5 



M /pc 3 for the environment of the early solar system, and iBrasser et al.l ( 20061 ) requires central 
densities higher than 10 3 M /pc 3 to reproduce Sedna's orbit. In terms of the mean density of the 
material the Sun would have interacted with in the cluster environment, the Sun would have had 
to have encountered a mean de nsity greater than 1 3 an d less than ~ 10 4 M /pc 3 to be consistent 
with our survey observations. iGutermuth et al.l (120051 ) map the volume density of three young 



embedded cluster regions (GGD 12-15, IRAS 20050+2720, NGC 7129). The peak densities of these 
regions were on the order of ~ 10 5 M /pc 3 . For GGD 12-15 and IRAS 20050+2720, 72% and 91% 
of the member stars reside in locations with densities upwards of 10 4 M /pc 3 , and are unlikely to 
produce the observed Sedna population. For NGC 7129, less than 24% of th e stars in the core of 
the cluster experience densities greater than 10 4 M /pc 3 . lLada &: Ladal (|1995l ) estimate the central 



stellar density of the 0.1 pc central regions of IC 348, NGC 2024, and Trapezium clusters to range 
from ~ 10 3 — 10 4 Mq/pc 3 3 at the minimum central density required to form Sedna's orbit. These 
environments and NGC 7129 could produce the observed Sedna population. 



5.4. Population Estimate 

Now that we have found the 10 M /pc 3 cluster population is the only cluster environment 
capable of emplacing Sedna on its orbit, we can place constraints on the size of the produced 
population. To estimate the size of the Sedna population, we use the value of a measured by Fraser 
et al (2010) for the hot and cold populations of the Kuiper belt (a=0.35 and 0.82 respectively) as 
limits for our brightness distribution. For each given value of a, absolute magnitudes are randomly 
assigned to our survey simulator created 3,000,000 Sednas 50,000 times, for every value of Nh<i.q- 
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A single instance of the brightness distribution can be thought of as a separate survey. For each 
Njki.6 tested, the number of synthetic "surveys" in which, like the real survey, one object on 
a Sedna-like body is detected are tallied. Valid detections are only those in which the object is 
located on the same CCD and in all 4 field observations. We do not require that the object have 
Sedna's absolute magnitude (H=1.6), only that the apparent magnitude of the object is above the 
SExtractor calculated limiting magnitudes of all 4 frames the object is "discovered" on. 

The best-fit values for the number of objects brighter than or equal to Sedna with 95% errors 
are 393 1 264 and 74+ 47 for the hot and cold brightness distributions respectively. The lower 
and upper 95% confidence levels limits reported are one-sided statistics found by computing the 
interval over which the integrated probability distribution 0.95 respectively of the total area. The 
survey simulator assumes all simulated Sednas that land on our images and are above the image 
limiting magnitude would be detected in the survey. The effect of a less than 100% survey detection 
efficiency is non-negligible. The reported size estimates represent a lower-bound on the size of the 
Sedna population. Assuming a uniform detection efficiency which drops to zero at the image 
limiting magnitude, the best-fit value and 95% limits for Nh<is is scaled by the inverse of the 
survey efficiency. For our nominal detection efficiency of 0.66, the best-fit values for the number of 
objects brighter than or equal to Sedna are 595l 400 and H2l 71 respectively for the hot and cold 
brightness distributions. Figure [9] plots the fraction of simulated surveys that produced a single 
Sedna detection as a function of Nh<i.g the 10 4 M Q /pc 3 cluster environment for our nominal survey 
detection efficiency. 

For the 10 4 M /pc 3 cluster environment, the range is quite large but there could be on the 
order of hundreds to thousands of planetoids brighter than Sedna present beyond the Kuip er belt. 



For c omparison, the total number of Sedna-sized or larger bodies in the Kuiper belt is ~5-8 (JBrown 



2008); there may be an order of magnitude or two more mass residing in the Sedna region than 
exists in the present Kuiper belt. The expected number of objects with H<1.6 varies significantly 
with the slope of the brightness distribution. Choosing a steeper power law for the brightness 
distribution decreases the likelihood of detecting only one Sedna because of the larger number of 
bright objects populating detectable orbits and decreases the best-fit number of objects brighter 
than Sedna. Selecting a smaller value of a, a shallower brightness distribution, increases the 
likelihood of detecting only one object on Sedna's orbit by decreasing the number of synthetic 
surveys with multiple detections. 

We excluded orbits with semimajor axes greater than 3000 AU from our analysis. If we had 
included higher semimajor axis orbits, the best-fit number of Sedna-like bodies brighter than or 
equal to Sedna would increase due to the larger population of orbits being included, but our overall 
conclusions would not change. Those orbits that are contributing most to being Sedna detections 
are those with semimajor axes much smaller than 3000 AU. The 10 4 M /pc 3 cluster results are the 
most sensitive to this cut. We examined the single detections produced in the best-fit simulations 
for the range of a parameters. Orbits with semimajor axes less than 1000 AU contribute ~80- 
of the Sedna single Sedna detections. The number is even greater from the 10 6 and 10 5 
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M0/pc 3 clusters. We make a conservative perihelia cut at 50 AU to compare the cluster produced 
Sedna-like orbits to our observed Sedna population. Our simulations were rerun, including orbits 
with perihelia greater than 45 AU and our conclusions remain the same. The 10 6 and 10 5 M /pc 3 
cluster will be ruled out with a higher confidence level because of the increase in low perihelia orbits 
that should have been detected in our survey. For this analysis we used the SEXtractor computed 
limiting magnitudes to determine whether our simulated Sedna population was observable on our 
images. 



5.5. Comparison to Occultation Surveys 



Wang et al.l (|2009i ) place upper limits on the number of small bodies in the Sedna region due 
to the lack of occultations from distant solar system bodies in the T AOS survey. We ca n estimate 
whether our upper limits for the Sedna population are consistent with lWang et al.l ()2009l )'s reported 
upper limits on the number density of objects larger than 1 km. We find the fraction of our 10 4 
M Q /pc 3 cluster survey simulator created 3,000,000 Sednas that are within 3° of the ecliptic (TAOS's 
ecliptic latitude range) and at 100 and 1000 AU. Assuming no break in the size distribution, we 
extrapolate the number of bodies larger than 1 k m. The albedo distribution is uncerta in, Sedna's V 
albedo is measured to be between 0.16 and 0.30 (|Brownll2008l : IStansberry et al.l 120081 ) . but to give 
an extreme upper limit we chose an albedo for 0.04 and assume no break in the size-distribution 
in order to estimate the fraction of bodies that would be observable by TAOS. TAOS is sensitive 
to bodies brighter than H=19.1. 

Within 3° of the ecliptic, 0.05% of the 10 M /pc 3 cluster produced Seda population are 
located between 50-150 AU and 0.07% reside at 900-1100 AU. For the flat size distribution value 
(a=.35), we expect there to be no more than 780 Sednas/deg 2 on the ecliptic at 100 AU and 10 4 
Sednas/deg 2 at 1000 AU assuming a 66% detection efficiency Our 95% confidence level estimates 
for a flat size distribution are well below TAOS's ecliptic number density of 1 km or larger bodies 
at 100 (~10 Sednas/deg 2 ) and 1000 AU (~10 9 Sednas/deg 2 ) even without a break in the size 
distribution. The TAOS observations do not rule out a large Sedna population with thousands of 
Sedna-sized or larger bodies residing far from the Sun for a flat brightness distribution. For the 
steep (cold population) size distribution, a = .82, and a 66% magnitude detection efficiency, at 
100 AU we expect no more than 2.8xl0 10 objects/deg 2 , approximately two orders of magnitude 
larger than TAOS's 3-<r upper limit. We find that even our 95% lower limit at 100 AU is an order 
of magnitude larger than the TAOS limit. Our expected number density at 1000 AU at our 95% 
upper confidence level, on the other hand, is 3.70xlO n objects/deg 2 below TAOS's limit of ~10 12 
objects/deg 2 . 

The occultation results do not necessarily rule out a steep size distribution for the Sedna 
population. In t he Kuiper belt at sma l l sizes (^50-150 km ) the distribution is observ e d to break to a 



shallo wer slope (IBernstein et al. 



2004 



Fuentes et al.ll2009l : iFraser fc Kavelaars!l2009l ). iBrasser et al. 



(|2006l )'s model did not include gas in the solar nebula and therefore did not include the effects of 
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gas dynamics in their simulations. Sedna is ^1500 km in size and would not be effected by gas 
drag, but smaller sized objects would be. iBrasser et al.l (|2007l ) investigated the effect of gas drag 
on the size distribution of objects deposited into the Sedna region. They find a size sorting effect 
in the cluster-produced Sedna population. Bodies smaller than ~20-60 km would be circularized 
onto orbits beyond Jupiter and Saturn and not available to be scattered into the Sedna region. 
Far fewer small-sized objects would be deposited into the Sedna region. Our survey is sensitive to 
objects much larger than those that would be effected by gas drag or the break in the brightness 
distribution. The combination of a broken power-law size distribution and a size-sorting effect 
could reconcile the observations, causing very few small objects that TAOS would have been able 
to detect to be present in the Sedna region. 



5.6. Open Cluster Environments 



The majority of stars are birthed in embedded clusters, b ut 4-7% of s tars form in smaller loose 
conglomerations with little or no gas known as open cl usters (ILadal 12004! ) . Open clusters, like the 
Pleiades, have ages of a few tens to hundreds of Myrs (JLadal 120041 ) . Although embedded clusters 
are more prevalent, it is postulated that ^ 5 % of the embedded c lusters may dissipate into loosely 
bound open clusters (ILada &: Ladal 120031 ) . iKaib fc Quinnl (120081 ) are able to produce objects on 
Sedna-like orbits in various open cluster environments. Interactions between the planetesimals disks 
of the cluster members are not included in th eir simulations . Thei r numerical integrations produce 
similar wedge-l i ke orb ital distributions to the lBrasser et al.l (|2006l ) embedded clusters models, but 
Kaib fc Quinnl (|2008l ) find no relationship between the size of the birth cluster and the orbital 
distribution of Sednas. The open cluster integrations are nondeterministic with Sedna's orbit being 
produced in only 5 of their 16 cluster simulations of varying cluster size. For those integrations 
that do produce Sed na and other Sedna-like orbits, distributions similar to the 10 4 and 10 6 M Q /pc 3 
Brasser et al.l (J2006J) results are generated. This is not unsurprising since the dominant dynamics 
sculpting the Sedna region, stellar encounters, is the sa me in both environmen ts. Our analysis above 
of the embedded cluster distributions also applies to IKaib &: Quinnl (120081 ) open cluster orbital 
distributions. Those distributions where Sedna is at the end of a distribution Sedna-like orbits 
with many lower semimajor axes and lower perihelia orbits similar to the 10 5 and 10 6 M0/pc 3 3 
embedded clusters, are inconsistent with our observations. 



5.7. Implications for the Kuiper Belt 



Using the discov ery of 2008 KV42, with an orbit essentially perpendicular to the ecliptic, 
Gladman et al.l (120091 ) posit a metastable parent population with inclinations greater than ~50 AU 
with a in the hundreds of AU and q = 3545 AU. Such a population is produced in t he 10 5 and 10 6 



M^/ pc 3 cluster environments but not present in the 10 4 Mq/pc 3 embedded cluster ( Brasser et al 



20061 ) . I Gladman et al.l (|2009l ) suggest 2008 KV42 may have been a high inclination counterpart 
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to Sedna placed on a lower perihelia and semimajor axis that later diffused to its current orbit. 
Although for our analysis we removed such objects with perihelia less than 50 AU from our dis- 
tribution, adding those objects would only rule out the 10 5 and 10 6 M /pc 3 cluster environments 
to even higher confidence because many more low a and low q objects single detections would be 
produced than detections with similar orbits to Sedna. With this region devoid of particles in the 
10 Mq/pc 3 cluster integrations, this suggests that 2008 KV42 and Sedna are likely formed from 
two independent source populations. 



6. Latitude Distribution 



Figure [TUl plots the folded latitude distribution of all objects wi th a>30 deb i ased f or latitudinal 
coverage. We assume Poisson detection statistics (as computed bv lKraft et al.l (|199ll )). with error 
bars representing the Poissonian 68% confidence limit on the detected number of objects in each 
latitud e bin corrected for sky coverage. A noticeable spike occurs at ~12° from the ecliptic. [Brown 
(J2008I ) also finds these prominent peaks in the latitudinal distribution ~± 11° ecliptic latitude. 
Brown finds that this peaked dist ribution canno t be generated by a simple inclination distribu- 
tion of objects in random orbits. iBrownl (|2008l ) suggests that resonant orbits are likely able to 
explain these high latitude concentrations. Resonators trapped in the Koz ai resonance (su ch as 
Pluto) have their perihelia near their maximum excursion off the ecliptic (jMorbidellil 119971 ) and 
the highest d etection probab i lity o ut of the ecliptic plane. The plutinos come to perihelia away 
from Neptune (|Malhotralll996l . Il995l ) and are pref e rentia lly biased towards detection at certain lon- 
gitudes. Without dynamical classification IBrownl (|2008l ) could not verify the plutinos as the source 
of these peaked latitude distributions. 

With secure orbits for our detections we can address this issue. In order to classify which 
of the survey KBOs reside in mean motion resonances with Neptune, each KBO had 13 clones 
integrated for 10 MYrs. One clone represents the best-fit orbit, and the rest are taken from a 
self-consistent spread of orbits covering the 3-<r uncertainty of the KBOs best-fit orbital solution 
computed from the covariance matrix of orbital elements obtained f rom AstDycl on 2009 Dec ember 
1. These objects were integrated using the n-body c ode SyMBA (ILevisqn &: Duncan! Il994l ) using 
the integrator swift_rmvs3 based on the mapping bv lWisdom &: Holmanl (|199ll ). The KBO clones 
were treated as massless particles. The four giant planets were included and their initial conditions 
were taken from JPL HORIZONS. The mass, position, and velocity of the terrestrial planets were 
combined with the Sun. The integration proceeded backwards in time with 40-day time steps from 
epoch JD 2455200. After 10 MYrs, the clones were examined for one or more librating resonant 
angles and as well as librating arguments of perihelion in order to identify Kozai resonators. We 
identified objects (listed in Table [2]) as resonant if all the clones lie in the resonance at the end of 
the integrations. 



http : //hamilton . dm . unipi . it/astdys/ 
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The latitudinal distribution of detected plutinos found in the survey is plotted in Figure [TTJ 
Six plutinos were detected in our survey, only two reside at ecliptic latitudes less than 10°. The 
remaining four plutinos compose the majority of the 12° latitude spike. Of these four plutinos 
(Huya, 2007 RT15, 2002 VE95 and 2008 SO2006), two objects are Kozai resonators, Huya and 2007 
RT15; the other three have perihelia off the ecliptic having possibly experienced temporary Kozai 
interactions. The remaining non-plutino distribution still exhibited a peak in the distribution at 12° 
including two members of the Hau mea collisonal fami l y. At least 7% of our detec t ions are fragments 
of the Huamea collisional fa mily ([Brown et al.l 120071 ; iRagozzine fc Brownl 120071 ; ISchaller fc Brown 
20081 : ISnodgrass et al.l |2010| ) . The identifier of the Haumea f amily is the characteristic deep near- 
infra red pure water ice absorption features on their surfaces ([Brown et al.l 120071 ; ISchaller fc Brown 
20081 ) . The water ice -rich bodies are thought to all have anomalously high albedos, like family 
member 2002 TX300 ([Elliot et al.ll2010l ). extremely biasing our survey toward detection of Haumea 
family members. Any clustering in the Haumea family members will severely bias our latitude 
distribution. Removing the spectroscopicly confirmed family members from our survey, the non- 
plutino distribution is not peaked as shown in Figure [TTJ 



Brownl (J2008I ) and this work are the only two wide-field surveys to probe significantly beyond 



the ecliptic. In order to test whether the plutino population observed by ecliptic surveys is represen- 
tative of the entire plutino population, we co mpare our observed p l utino latitude distrib ution to the 



CFEPS plutino model. The CFEPS survey (JKavelaars et al. 



2009 



Gladmanet al.ll2010i ) orbital and 



brightness distribution is based on the sample of plutinos detected in observations covering ecliptic 
latitudes less than 2°. None of their detections are Kozai librators, thus only representing the 
non-Kozai plutino population. CFEPS is sensitive to an absolute magnitude range of H g /~ 6-10.5, 
fainter than the sources we are able to detect in our survey. In order to compare their model to our 
observed latitude distribution, we must extend the distribution to larger objects where the CFEPS 
survey does not measure directly and where the slope of their measured brightness distribution may 
not be applicable to the larger bodies that w e detec t. The H distribution is measured in g' and we 
observe in a broadband R filter. iFraser et al.l ([20081 ) find an average KBO value of < g'-R >= 0.95, 
and we apply this as our constant offset to the g' magnitudes. We create a latitude distribution 
by shuffling the absolute distribution of the 10 5 model plutinos with H 9 >10.5 ~10 5 times. We 
tally the latitudes of all plutinos for all runs with magnitudes brighter than R=22 that lie within 
our survey sky coverage in a folded latitude histogram binned in 2° latitude bins. To estimate the 
expected number of plutinos in the 12° bin, we scale CFEPS model latitude distribution to the 
value of our folded latitude distribution at 2°, the lowest latitude binned plutino detection in our 
survey. Assuming Poisson errors and using the quadrature of the fractional errors. 6.7 _ 6 j (68% 
confidence level) times as many plutinos reside in 11-13° from the ecliptic than are predicted by 
the non-Kozai plutino CFEPS model. Although the range is quite high, our latitude distribution 
suggests that the plutino population in particular the Kozai population has been underestimated 
and may be much larger than previous KBO surveys have reported. 
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7. Number of Bright Objects 

Our survey probes the bright-end of the KBO size distribution. Assuming a uniform latitude 
distribution we can crudely estimate the number of large observable KBOs. Using our nominal 
survey efficiency function from Section f3. 3, 21 and the effective area covered we compute the expected 
number of bright KBOs (a>30) as a function of magnitude. We neglect the effects of masked CCD 
regions and other geometric effects. Our sky coverage drops significantly above latitudes of ±30° 
and therefore we only focus on detections and sky covered within ±30° of the ecliptic. Figure [12] 
plots the cumulative number of expected bright KBOs as a function of magnitude compared to 
known multiopposition KBOs. 



8. Conclusions 

Surveying ~12000 deg 2 within ±30° of the ecliptic to ~21.5 in R magnitude, we have searched 
for additional members of the Sedna population. Based on the 52 KBOs and Centaurs detected in 
our survey we conclude: 

• We detected only one object on a Sedna- like orbit, Sedna, despite a sensitivity to motions 
of bodies out to ~1000 AU. With one detection, we cannot differentiate between the various 
proposed formation mechanisms proposed to emplace Sedna on its orbit. 

• For the embedded cluster Sedna formation model, we reject the 10 5 and 10 6 M Q /pc 3 cluster 
environment-produced populations as consistent with our redetection of Sedna. We find the 
10 4 M /pc 3 cluster environment consistent with our observations, with a best-fit population of 
Nh<i.6= 595^ 400 for the hot population and 112l 71 for the cold population size distributions 
assuming our nominal detection efficiency of 66%. 

• The plutino population has a peaked distribution at ~±12° ecliptic latitude, likely due to 
Kozai resonators and current estimates of the size of the plutino population from on-ecliptic 
surveys insensitive to these high latitude plutinos likely underestimate the size of the true 
population 
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Table 1:: Summary of Field Positions 



Pointing 


R.A. 

(J2000) 


Dec. 

(J2000) 




Night 1 






Night 2 






MJD obsl 


MJD obs 2 


mag. limit 


MJD obs 1 


MJD obs 2 


mag limit 


1 


13 53 27.599 


-18 15 59.80 


54228.327 


54228.333 


20.9 


54229.189 


54229.237 


21.4 


2 


13 55 33.959 


-18 15 59.80 


54228.330 


54228.336 


21.1 


54229.192 


54229.240 


21.5 


3 


13 22 04.440 


-13 43 28.20 


54228.313 


54228.320 


20.6 


54229.163 


54229.209 


20.9 


4 


13 38 45.960 


-13 43 28.60 


54228.260 


54228.266 


21.4 


54229.169 


54229.216 


21.3 


5 


13 40 49.440 


-13 43 27.50 


54228.263 


54228.269 


21.4 


54229.172 


54229.219 


21.4 


6 


13 55 27.841 


-13 43 27.80 


54228.205 


54228.211 


21.2 


54229.176 


54229.222 


21.3 


7 


12 51 16.560 


-09 10 55.60 


54228.273 


54228.279 


21.0 


54229.183 


54229.230 


21.4 


8 


12 53 18.240 


-09 10 55.90 


54228.276 


54228.283 


21.1 


54229.186 


54229.233 


21.4 


9 


13 07 40.440 


-09 10 55.90 


54228.219 


54228.226 


21.3 


54229.196 


54229.243 


21.5 



Table 1:: Full table can be found in the online version. The 
center coordinates for all pointings searched for KBOs and 
used in the analysis presented in this paper. The table in- 
cludes pointing number, the right ascension and declination 
of the mosaic center CCD (B15), MJD dates of all four ob- 
servations of the field and limiting magnitudes for each night 
the field was observed 



to 



27 



designation 


a 


c 


i 


R 


oppositions 


night 1 


night 2 


H 


MMR 




(AU) 




(deg) 


(AU) 




avg mag 


avg mag 






(26181) 1996 GQ21 


93.01 


0.588 


13.4 


40.87 


11 


21.0 


20.6 


5.2 


11:2 


(26308) 1998 SM165 


47.99 


0.375 


13.5 


36.99 


12 


21.4 


21.4 


5.8 


2:1 Kozai 


(19521) 1998 WH24 


45.74 


0.103 


12.0 


41.84 


11 


21.9 


21.9 


4.8 




(40314) 1999 KR16 


48.83 


0.306 


24.8 


36.27 


6 


20.6 


20.7 


5.8 




(38628) 2000 EB173 


39.44 


0.277 


15.5 


28.93 


7 


19.3 


19.4 


4.7 


3:2 Kozai 


(47932) 2000 GN171 


39.36 


0.281 


10.8 


28.34 


9 


20.3 


20.4 


6.0 


3:2 


(20000) 2000 WR106 


42.85 


0.056 


17.2 


43.39 


13 


19.9 


19.7 


3.6 




(82075) 2000 YW134 


57.61 


0.287 


19.8 


43.81 


6 


21.0 


21.1 


4.9 


8:3 


(83982) 2002 G09 


19.45 


0.277 


12.8 


15.00 


6 


21.1 


21.1 


9.1 




(50000) 2002 LM60 


43.47 


0.039 


8.0 


43.27 


16 


19.3 


19.4 


2.5 




(55636) 2002 TX300 


43.46 


0.126 


25.8 


41.29 


12 


19.8 


19.7 


3.3 




(55638) 2002 VE95 


39.37 


0.290 


16.3 


28.24 


10 


20.1 


20.0 


5.3 


3:2 


(119979) 2002 WC19 


47.80 


0.260 


9.2 


42.95 


7 


21.4 


21.4 


5.0 


2:1 


(174567) 2003 MW12 


45.87 


0.144 


21.5 


47.95 


12 


20.6 


20.5 


3.6 




(120178) 2003 OP32 


43.45 


0.108 


27.1 


41.36 


6 


19.9 


19.9 


4.1 




(120181) 2003 UR292 


32.49 


0.176 


2.7 


26.87 


6 


21.4 


21.5 


7.0 




2003 UZ117 


44.29 


0.133 


27.4 


39.46 


6 


22.0 


21.7 


5.3 




(90377) 2003 VB12 


510.00 


0.850 


11.9 


88.31 


8 


21.1 


21.0 


1.6 




(136204) 2003 WL7 


20.17 


0.259 


11.2 


15.21 


6 


20.6 


20.6 


8.7 




(175113) 2004 PF115 


39.18 


0.062 


13.4 


41.34 


6 


20.5 


20.3 


4.7 




2004 PG115 


92.08 


0.605 


16.3 


36.65 


5 


20.8 


20.9 


5.0 




(120347) 2004 SB60 


42.27 


0.105 


23.9 


43.89 


10 


20.4 


20.5 


4.2 




2005 CB79 


43.15 


0.140 


28.7 


40.16 


6 


20.4 


20.3 


5.0 




(145451) 2005 RM43 


91.37 


0.616 


28.8 


35.19 


7 


19.9 


19.9 


4.4 




(145452) 2005 RN43 


41.77 


0.028 


19.2 


40.72 


13 


20.0 


20.0 


3.9 




(145480) 2005 TB190 


76.58 


0.397 


26.4 


46.45 


7 


20.9 


20.8 


4.7 




2006 SX368 


22.28 


0.463 


36.3 


12.44 


4 


20.3 


20.3 


9.5 




2007 JF43 


39.41 


0.185 


15.1 


39.45 


4 


20.9 


20.8 


5.2 


3:2 


2007 JF45 


44.69 


0.147 


10.6 


38.12 


Id 


21.5 


21.4 


6.0 




2007 JJ43 


48.22 


0.166 


12.0 


41.96 


3 


20.8 


20.7 


4.9 




2007 JK43 


46.35 


0.492 


44.9 


23.93 


3 


20.8 


21.1 


7.6 




2007 NC7 


34.39 


0.507 


6.3 


20.37 


3 


21.4 


21.6 


8.6 




2007 OC10 


50.09 


0.292 


21.7 


35.48 


3 


20.8 


20.8 


5.7 




(225088) 2007 OR10 


67.34 


0.500 


30.7 


85.37 


7 


21.5 


21.4 


1.9 




2007 PS45 


43.75 


0.090 


18.9 


39.80 


Id 


21.5 


21.1 


5.6 




2007 RG283 


19.98 


0.233 


28.8 


18.70 


3 


21.5 


21.0 


8.8 




2007 RH283 


15.96 


0.339 


21.4 


17.48 


8 


21.4 


21.2 


8.4 





Continued on Next Page. . . 
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designation 


a 


c 


i 


R 


oppositions 


night 1 


night 2 


H 


MMR 




(AU) 




(deg) 


(AU) 




avg mag 


avg mag 






2007 RT15 


39.61 


0.234 


12.9 


30.90 


3 


21.6 


21.3 


6.9 


3:2 


2007 RW10 


30.40 


0.303 


36.0 


26.24 


7 


21.3 


21.1 


6.5 




(229762) 2007 UK126 


73.52 


0.488 


23.4 


45.96 


9 


20.4 


20.3 


3.4 




2007 XV50 


46.02 


0.073 


22.9 


46.19 


3 


21.2 


21.3 


5.0 




2008 AP129 


41.66 


0.138 


27.4 


37.39 


5 


20.6 


20.7 


5.3 




2008 CS190 


42.08 


0.153 


16.0 


36.17 


2 


21.6 


21.6 


6.4 


5:3 


2008 CT190 


52.47 


0.339 


38.9 


34.77 


2 


21.0 


21.4 


5.5 


7:3 


2008 LP17 


88.04 


0.660 


14.1 


30.26 


2 


21.0 


20.9 


6.6 




2008 NW4 


45.58 


0.203 


23.1 


36.92 


2 


21.2 


21.0 


6.0 




2008 OG19 


67.37 


0.428 


13.1 


38.74 


2 


21.6 


21.3 


4.9 




2008 QB43 


43.36 


0.219 


26.3 


38.79 


3 


21.6 


21.4 


5.6 




2008 QY40 


63.09 


0.418 


25.1 


38.11 


2 


20.9 


20.9 


5.3 




2008 S0266 


39.64 


0.247 


18.8 


31.58 


2 


21.5 


21.4 


6.9 


3:2 


2008 SP266 


41.21 


0.124 


19.5 


36.18 


2 


21.2 


21.2 


5.7 




2008 ST291 


106.00 


0.607 


20.7 


56.68 


2 


21.8 


21.3 


4.4 





Table 2:: Orbital elements reported by the Minor Planet Cen- 
ter of Centaurs and KBOs detected in the Palomar survey: 
semimajor axis (a) , eccentricity (e) , inclination (i), barycen- 
tric distance (R), oppositions observed (in years excepted 
where days noted by d), nightly discovery magnitudes, ab- 
solute magnitude (H), and mean motion resonance (MMR) 
if applicable. All 2007 and 2008 classified objects were new 
discoveries found by this work 



cluster central density 






a 






{M Q /pc 3 ) 


0.2 


0.35 


0.4 


0.6 


0.82 


10 4 


60 


54 


48 


40 


47 


10 5 


99 


98 


99 


98 


97 


10 6 


100 


100 


100 


100 


100 



Table 3: 3D KS test results for the lBrasser et al.l (J2006J) cluster produced single detections compared 
to Sedna's orbit. We report the confidence level at which we can reject the two distributions as 
drawn from the same parent population. 
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Fig. 1. — Scale drawing of the focal plane of the QUEST camera, depicting the layout of the 112 
CCDs. The gap between CCDs in the north-south direction is ~1.2' and the spacing between 
adjacent fingers along the east-west direction is ~25'. 
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Fig. 2. — Sky coverage of the Palomar survey plotted on the J2000 sky. The observed fields are 
plotted to scale. The plane of the Milky Way is denoted as a dashed line, and the ecliptic is denoted 
as a solid line. Holes are due to galactic plane avoidance, bad weather, forest fires, and hardware 
malfunctions. 
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Fig. 3. — Survey efficiency based on the previously known multiopposition KBO population. Solid 
line plots the best-fit efficiency function. Diamonds plot the binned detection efficiency in 0.5 
magnitude bins with one-c Poissonian error bars. The dashed line with triangles is the number 
of previously known multiopposition KBOs (a i 30 AU) in each magnitude bin with predicted 
positions located on all 4 survey images. 
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Fig. 4. — Effect of geometric losses on the sky coverage of the survey. The fraction of simulated 
orbits found on all four CCDs as a function of latitude binned in two degree bins. Main effects 
are due to KBOs that are not located on all field observations and move off the CCD or objects 
positioned in the gaps between the CCDs. 
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Fig. 5. — Efficiency of our orbit-fit detection filter. Fraction of synthetic orbits with barycentric 
distances between 15-1000 AU that successfully are identified as outer solar system bodies found 
on all four survey images versus semimajor axis binned in 100 AU bins 
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Fig. 6. — Eccentricity vs. Semimajor axis and Inclination vs. Semimajor axis of multiopposition 
obje cts found in the Palomar survey . Sedna has been excluded for better resolution. One-cr errors 
from [Bernstein &: Khushalanil ([2000) orbit fit are plotted. The error bars are typically smaller than 
the size of the symbol. 
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Fig. 7. — Inclination vs. barycentric distance at discov ery for objects detected in the Palomar 
survey. One-u errors from iBernstein fc Khushalanil (I2000J ) orbit fit are plotted. The error bars are 
typically smaller than the size of the symbol. 
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Fig. 8. — Plot of perihelion (q) vs. semimajor axis (a) a nd plot of incl i nation (i) vs. semimajor 
axis (a) for Sedna-like bodies produced at the end of the iBrasser et al.l (J2006I ) embedded cluster 
simulations used in this work. We limit the population to orbits q > 50 and a < 3,000 AU. The 
diamond (red in the online version) denotes Sedna's orbit. 
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Fig. 9. — Results of the 10 4 M Q /pc 3 cluster analysis for the a=0.35 (hot) and a=0.82 (cold) Kuiper 
belt population size distributions- Fraction of synthetic surveys with one detectable Sedna-like body 
as a function of the number of bodies bigger and brighter than Sedna assuming our nominal 66% 
detection efficiency. 
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Fig. 10. — The folded latitudinal distribution of objects with semimajor axis > 30 AU found in 
this work. The lower dashed line with diamonds shows the number of actual KBO detections in 
two-degree bins. The dashed line shows the fractional ecliptic latitude completeness. The solid 
line shows the expected number of KBOs brighter than 21.3 corrected for sky coverage with one-cr 
Poisson error bars computed for the unfolded distribution added in quadrature. 
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Fig. 11. — The folded latitudinal distribution two-degree bins of survey of multi-opposition plutinos 
(top) and non-plutinos excluding confirmed Haumea family members (bottom) brighter than 21.3 
corrected for sky coverage with one-<7 Poisson error bars computed for the unfolded distribution 
added in quadrature. Filled circles shows the number of actual KBO detections in two-degree bins. 
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Fig. 12. — Cumulative number of expected KBOs within ± 30° of the ecliptic assuming a flat 
latitude distribution (open diamonds) with 2-cr Poisson error bars. Cumulative number of known 
multiopposition KBOs (a> 30)within ±30° ecliptic latitude (filled circles) 



